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We present measurements of the branching fractions for the Cabbibo suppressed decays D + —> 7r + 7r° 
and D + — > K + tv° based on a data sample corresponding to an integrated luminosity of 124.3 
fb _1 . The data were taken with the BABAR, detector at the PEP-II B Factory operating on and 
near the Y(4S) resonance. We find B{D+ -> tt+tt ) = (1.25 ± 0.10 ± 0.09 ± 0.04) x 10~ 3 and 
B{D+ -» K+tx°) = (2.52 ± 0.47 ± 0.25 ± 0.08) x 10~ 4 , where the first uncertainty is statistical, 
the second systematic and the last error is due to the uncertainties in the absolute branching 
fraction scale for D + mesons. This represents the first observation of the doubly Cabibbo-suppressed 
D + — > K + tt° decay mode and a new measurement of the D + — > 7r + 7r° branching fraction. 

PACS numbers: 13.25.Ft, 11.30.Hv, 13.30.Eg 



Measurements of rare hadronic D + decays provide in- 
sight into SU(3) flavor symmetry, QCD dynamics, and 
weak flavor mixing Q] . Studies of these decays are useful 
for D°D° mixing analyses, which benefit from improved 
measurements of D + — > 7r + 7r° and D + — > K + K° branch- 
ing fractions in order to understand the size of the SU (3)- 
violating effects in D meson decays. In addition, doubly 
Cabibbo-suppressed decays such as D° — ► K + ir~ com- 
plicate measurements of flavor oscillations in hadronic 
D° decays. Knowledge of the S , C/(3)-related channels re- 
ported here can lead to a better understanding of this 
background. Previous analyses of these D + decays were 
reported by MARK III, CLEO, and FOCUS jllEll- 

This analysis is based on data recorded with the BABAR 
detector at the PEP-II asymmetric-energy e + e _ storage 
ring at the Stanford Linear Accelerator Center. The 
data sample corresponds to an integrated luminosity 
of 124.3 fb _1 recorded at center-of-mass (CM) energies 
y/s = 10.58 GeV and 10.54 GeV and includes approxi- 
mately 167 x 10 6 e + e _ — > cc events. 

The BABAR detector is described in detail elsewhere 5] . 
Charged particle momenta are measured with a 5-layer 
double-sided silicon vertex tracker (SVT) and a 40-layer 
drift chamber (DCH), both inside the 1.5 T magnetic 
field of a superconducting solenoid. A calorimeter (EMC) 
consisting of 6580 CsI(Tl) crystals measures electromag- 
netic energy. Charged hadron identification is provided 
by measurements of the rate of ionization energy loss, 
dE/dx, in the tracking system and of the Chercnkov 
angle obtained from a ring-imaging Cherenkov detector 
(DIRC). The instrumented flux return of the magnet al- 
lows discrimination of muons from pions. 

We use a Monte Carlo simulation of the BABAR detec- 
tor based on GEANT 4 Q to validate the analysis and to 
determine the reconstruction efficiencies for the two sig- 
nal modes D + — > tt + ti° and D + — ► K + ir°, as well as for 
the decay D + — > K ~7r + 7r + , which is used as a reference 



to normalize our results. Simulated events are generated 
with the Pythia event generator [tJ. 

We reconstruct D + meson candidates in the signal 
modes by combining a charged track, identified either as 
a pion or kaon, with a reconstructed 7r° candidate. Until 
better knowledge of the overall e + e _ — > D + X produc- 
tion rate is obtained, any measurement of absolute D + 
branching fractions will be limited by the uncertainty in 
the number of D + mesons in the data sample. We avoid 
this uncertainty by measuring our signal modes relative 
to the high statistics, well-measured D + — ► K~n + ir + 
decay mode. 

In order to reduce the large amount of combinatorial 
background in the D + signal modes, we include only D + 
mesons that originate from D* + — ► D + ir° decays. To 
minimize systematic uncertainties in the reconstruction 
of the low momentum ir° from the D* + decay, this is 
done for both the signal and reference channels 

Only events with at least three charged tracks are se- 
lected for this analysis. Charged tracks are required to 
have a distance of closest approach to the interaction 
point in the plane transverse to the beam axis of less 
than 1.5 cm, a distance of closest approach along the 
beam direction of less than 10 cm, a minimum transverse 
momentum of 100 MeV/c, and at least 12 DCH hits. All 
candidate tracks in the reconstructed decay chains must 
satisfy a set of pion or kaon identification criteria based 
on the response of the DIRC and the dE/dx measure- 
ments in the tracking system. 

A pair of energy clusters in the EMC, which are iso- 
lated from any charged tracks and have the expected lat- 
eral shower shape for photons, is considered a ir° candi- 
date if both clusters exceed 30 MeV, and the associated 
invariant mass of the pair is between 0.115 GeV/c 2 and 
0.150 GeV/c 2 . The energy of the 7r° candidate in the 
laboratory frame is required to be greater than 0.2 GeV. 



Unless explicitly stated, charge conjugate reactions are implicitly 
included throughout this paper. 



We accept a D + candidate if its invariant mass falls 
between 1.7 GeV/ c 2 and 2.0 GeV/ c 2 . In addition we re- 
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quire that the cosine of the helicity angle, &h 2 , which is 
uniformly distributed for signal events but peaks at ±1 
for background, satisfies —0.9 < cos(0/,) < 0.8 for the 
D+ tt+tt mode and -0.9 < cos(6 h ) < 0.7 for the 
D+ K+tt° mode. 

Two positively charged pion tracks and a negatively 
charged kaon track are combined in a vertex fit to form a 
D + — > K~ir + ir + candidate. We require the chi-squared 
probability of the vertex fit to be V > 0.001. The candi- 
date is accepted if the invariant mass of the D + lies be- 
tween 1.75 GeV/c 2 and 1.95 GeV/c 2 . The smaller range 
compared with the signal modes reflects the better reso- 
lution for this decay mode, which has only charged tracks 
in the final state. 

D + candidates are combined with a reconstructed 
tt° to select D* + — > D + tt° decays. An additional re- 
quirement on the center of mass momentum, pcm < 
0.45GeV/c, is applied to tt° candidates used in the D* + 
reconstruction. Only D* + candidates with a mass differ- 
ence Am = m D *+ — m D + less than 0.155 GeV/c 2 are ac- 
cepted for this analysis. A requirement on the normalized 
momentum of the D* + meson 3 , xd» > 0.6, correspond- 
ing to a fl* center of mass momentum greater than 2.9 
GeV/c, eliminates backgrounds from B meson decays and 
further reduces the combinatorial background. If more 
than one D* + candidate is reconstructed in an event, we 
choose the one with the larger xd* value. 

With these requirements applied to Monte Carlo 
events, we obtain reconstruction efficiencies of 7.8% for 
the D+ tt+tt mode, 5.9% for the D+ -> K+n° mode, 
and 8.5% for the D+ -> K--K+TT+ mode. 

We extract the signal yield for each of the three decay 
modes from the invariant mass distribution of the D + 
candidates with unbinned extended maximum likelihood 
fits. 

For the reference mode, the remaining background in 
the K~ir + ir + invariant mass distribution is described by 
a first order polynomial; the signal lineshape is modeled 
by a double Gaussian function. In order to accommodate 
possible differences in the D + momentum distribution in 
data and Monte Carlo events, a weight function describ- 
ing the relative change in reconstruction efficiency as a 
function of the D + momentum is included in the like- 
lihood function. Similar corrections are applied to the 
signal mode fits described below. A second weight func- 
tion is used in the D + — ► K~tt + it + fit to correct for 
potential differences between the K~it + tc + Dalitz plot 



2 We define 9^ as the angle between the direction of the charged 
daughter particle of the D + decay and the direction of the D*+ 
meson evaluated in the D + rest frame. 

3 The normalized momentum i^. is defined as xp* = 

Pd* 1 I \j s I^~ m \>- where p£i f is the momentum of the D* me- 
son in the CM frame and s is the square of the energy of the 
initial e + e — system. 



structure in data and simulated events. 

Not all of the D + candidates in our sample originate 
from D* + — > D + ir° decays. Some D + mesons from other 
sources can combine with a random 7r° in the event to 
pass the Am requirement. This background is uniformly 
distributed in Am; in the invariant mass distribution it 
peaks at the D + mass. While it might be expected that 
this effect cancels when we calculate branching fraction 
ratios between signal and reference modes, Monte Carlo 
studies indicate that this is not the case. A correction, 
extracted from data, must be applied to compensate for 
the difference in relative efficiencies between this peaking 
D + background in signal and reference modes, caused 
mostly by the helicity angle requirement. We use the 
Am sideband, shown in Fig. ^ to determine the cor- 
rected D + yield. The Am signal region is defined as a 
2er window around the nominal D* + — D + mass differ- 
ence; the sideband extends from 5cr above the nominal 
value to 0.155 GeV/c 2 . The signal contains D + decays 
from all sources, while D + mesons present in the Am 
sideband come from sources other than D* + — > D + w° de- 
cays. The invariant mass distributions for the signal and 
sideband regions are fitted simultaneously with identical 
signal shapes. The yield in the sideband is constrained 
to be greater than or equal to zero to avoid an unphys- 
ical enhancement of the signal yield should the peaking 
D + background fluctuate low. We scale the D + yield 
from the sideband by the ratio of the integrals under the 
combinatorial background curve in the signal and side- 
band regions and subtract this value from the yield in the 
signal region to extract the net yield of D* + — > D + ir°, 
D + — ► K~7r + ir + (tt + tt°, K + 7T°) decays corrected for the 
peaking background. 

For the D + — ► tt + tt and D + — > K + ir° signal modes, 
the double Gaussian signal function is replaced by a 
bifurcated Gaussian function, which gives a better de- 
scription of the increased width of the signal toward 
lower masses, caused by radiative losses associated with 
7T° reconstruction. An exponential function is added to 
the linear background parameterization to model back- 
grounds from misreconstructed decays such as D + — » 
K° S K+ with K% -> 7r°7r° or D° -> K%<k° with K% -> 
7r + 7r~ , where two of the three decay products are used to 
reconstruct a signal mode candidate. This background 
contribution peaks at low mass, but has a long tail that 
extends into the D + signal region. 

The D + — * K + n Q mode has an additional background. 
Decays of Df mesons to K + K a s final states with the 
Kg decaying to two neutral pions can cause an excess in 
the invariant D + mass distribution if they are mistak- 
enly reconstructed as D + — > K + ir°, and are combined 
with a random tt° to mimic a D* + signal. Based on a 
Monte Carlo study of Df — > K + K s events we model 
this additional background component with a Gaussian 
function centered at 1.807 GeV/c 2 and with a width of 
0.028 GeV/c 2 . 
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FIG. 1: Am distribution for simulated D* + -> D + ir°, D + -> 
if _ 7r + 7r + events. The vertical lines mark the signal and side- 
band regions defined in the text. The solid line shows the 
result of the fit; the background is given by the dash-dotted 
line. 



Signal and background shapes used in the fits are de- 
rived from Monte Carlo events. We minimize systematic 
uncertainties due to differences between data and the 
simulation by allowing most parameters to vary in the 
data fits. The only exceptions are in the D + — > A + tt° 
mode where the expected yield is too small to determine 
the signal shape parameters directly from data. Instead, 
we use the parameters found in the D + — > tt + tt data 
fit with the widths reduced by 5%. This correction was 
obtained in a Monte Carlo study of D + -> K + ir° and 
D + — > 7r + 7r events. The second exception is the shape 
of the Df — > K + Kg background which is constrained to 
values obtained from Monte Carlo simulations. 

A Monte Carlo sample corresponding to an integrated 
luminosity of approximately 80 fb _1 is used to validate 
the fit procedure. The branching fraction ratio of D + — > 
tt+tt to D + — > A~tt+tt+ decays in the simulation is 
2.8 x 10~ 2 . With the yields, Nf it , extracted from the 
D + — > tt+tt and D + — > A~tt+tt+ fits and the previously 
determined reconstruction efficiencies e, we obtain 

B(D+ -» TT+TT ) _ N flt (TT+7r°)-e K -^ + 

B(D+ -> K~TT+7T+) N flt {K-TT+TT+)-e^ +7T U 

= (2.7 ±0.1) x 10~ 2 (2) 

in good agreement with the expected value. As an ex- 
ample, the D + — > 7r + 7T° fit for simulated data and the 
different fit components are shown in Fig. [3 We repeat 
this study for the D + — ► K + tt° mode with the branching 
fraction ratio in the Monte Carlo sample set to 2.8 x 10 -3 
and find a value of (3.1 ± 0.6) x 10 -3 for the ratio of the 



FIG. 2: 7r + 7r u invariant mass distribution in the D + — > tv + tt° 
signal mode for simulated data. The background is modeled 
by a first order polynomial (dash-dotted) and an exponential 
(dotted), while the signal is fitted by a bifurcated Gaussian 
function. The dashed line shows the peaking D + background 
determined from the Am sideband fit. 



efficiency-corrected yields returned by the D+ — » K + 7r a 
and D+ -> K~-k+-k + fits. 

Figure shows the fit results for the full data sample. 
The signal yields are N flt {I<- tt+tt+) = 101380 ± 415 
for the reference mode, N ^(tt+tt ) = 1229 ± 98 for the 
D+ -> tt+tt signal mode and N flt (K+TT°) = 189.0 ± 
35.2 for the D + — > K + ir° signal mode. The errors are 
statistical only. 

The systematic errors in this analysis include uncer- 
tainties in the reconstruction efficiencies as well as the 
errors associated with the event yields returned by the 
maximum likelihood fit. Since we measure the branch- 
ing fraction for D + — > 7r + 7r° and D + — > K + ir° rela- 
tive to the D + — > K~tt + it + reference mode several sys- 
tematic uncertainties in the efficiencies cancel or are re- 
duced. Individual systematic contributions listed in the 
following paragraph apply to the measurement of the 
B(D+ -► tt+tt ) to B{D+ -> K-TI+TT+) ratio. Uncer- 
tainties of the B(D+ -> K+ir°) to B(D+ -> A'-tt+tt+) 
ratio are indicated in parentheses if they are different. 

The relative systematic error on the efficiencies in- 
cludes contributions of 1.9% from charged track recon- 
struction and vertexing, 0.9% (0.4%) from particle iden- 
tification, 3.2% due to uncertainties in tt° reconstruction, 
and 1.1% (1.2%) because of limited Monte Carlo statis- 
tics. Uncertainties in reconstruction efficiency of the slow 
tt° from D* + decays cancel in the branching fraction ra- 
tios. However, because of the 7r in the final state of 
the signal modes the D + mass resolution differs between 
signal and reference modes, which has an effect on the 
shape of the peak in the Am distribution. Based on a 
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study of simulated events we assign a 5% systematic er- 
ror due to this effect. The difference between the Am 
shapes in data and Monte Carlo events leads to an ad- 
ditional systematic uncertainty of 1.4% for the efficiency 
ratio. We vary the signal and background parameter- 
izations used in the fits of the D + invariant mass dis- 
tributions and derive a systematic uncertainty of 1.5% 
for the D + — ► 7r + 7r° to D + — ► K^n + n + efficiency ratio. 
The fixed signal shape used to determine the amount 
of peaking D + background in the D + — > K + tt° signal 
mode causes an additional systematic uncertainty that 
we estimate by varying the parameters of the Gaussian 
fit function within their errors and refitting the data. 
The quadrature sum of both effects gives a 2.7% system- 
atic uncertainty in the D + — > K + ir° to D + — ► K~n + ir + 
efficiency ratio. 

We add the individual contributions in quadrature and 
obtain a total systematic uncertainty for the D + — > tt + tt° 



to D+ 



K 7r + 7r + and D + 



K+ir° to D+ 



K tt+tt + efficiency ratios of 6.7% and 7.0%, respectively. 

The various weights, as well as the signal to sideband 
scale factor used in the maximum likelihood fits, give rise 
to systematic uncertainties in the D + yield. We study 
effects caused by the weight functions by varying their 
parameterizations within their errors, and find that this 
leads to a systematic error of 0.8% for the two signal 
modes and 1.2% for the D + — > A'~7t+tt+ reference mode. 

The scale factor used for the peaking D + background 
correction, extracted from data using the Am distribu- 
tion in the reference mode, is found to be (34.5 ± 1.6)%. 
The same scale factor is used for all three D + decay 
modes which results in an additional systematic error 
of 1.5%. For the D+ -> tt+tt and D+ -► K-tt+tt+ 
modes, we combine this in quadrature with the statisti- 
cal error of the scale factor to obtain a total systematic 
uncertainty in the amount of peaking D + background of 



4.9%. In the D+ -> K+tt° sample, the D+ yield in the 
sideband region fluctuates to a negative value; thus the 
nominal data fit constrains this component to zero. To 
estimate the uncertainty due to the peaking D + back- 
ground in this mode we repeat the fit without the con- 
straint and find a difference in signal yield between the 
constrained and the unconstrained fit of 9.5 events. We 
use this value as the systematic error on the background- 
corrected D + — > K + ir Q yield. 

An additional systematic uncertainty in the D + — > 
K +tt° analysis is due to the D+ — > K + Kg background 
component. Based on the measured branching fraction 
for this decay we expect 87 ± 58 D+ -> K+K° s back- 
ground events in our data sample. This is consistent with 
the fitted yield of 118 ± 53. When we vary the param- 
eterization of the Df — > K + Kg background function in 
the fit, the reconstructed D + — > K + ir° yield changes by 
4.5%, which is taken into account as a systematic error. 

With these systematic errors combined in quadrature, 
the yield in the reference mode is iV D +^^- 7r + 7r + = 
101380 ± 415 ± 1374. For the signal modes we find 
Nd+^+ttO = 1229 ± 98 ± 10 and N D+ ^ K +^o = 189.0 ± 
35.2 ± 12.8. 

Following Eq. we combine these measurements with 
the reconstruction efficiencies to obtain 

B{D+ -» 7T+7T ) 



B(D+ -> K-tt+it^ 



and 



B(D+ -> K+tt°) 
B(D+ -> K-K+-K+) 



= (1.33 ±0.11 ±0.09) x 10" 



= (2.68 ±0.50 ±0.26) x 10" 



Using B(D+ -> K-TT+7T+) = 0.094 ± 0.003 which is the 
weighted average of a recent CLEO-c result and the 
PDG value , we derive the branching fractions for the 
two signal modes 

B(D + -> tt+tt ) = (1.25 ± 0.10 ± 0.09 ± 0.04) x 10" 3 
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and 



B(D+ -» K 



(2.52 ± 0.47 ± 0.25 ± 0.08) x 10" 



where the last error is due to the experimental uncer- 
tainty in the D + — > K~it + it + branching fraction. Wc 
compute the significance S of the D+ — > K + ir° sig- 
nal as S = v / 2(ln£(N s ) - lnC(N s = 0)), where £(JV S ) 
is the maximum likelihood at the nominal fit yield, and 
C(N S = 0) is the value of the likelihood for N s = 0. We 
include systematic uncertainties by repeating this proce- 
dure while varying the fit parameters within their errors. 
The smallest signal significance obtained in this manner 
is 6.5 standard deviations. 

This represents the first observation of the doubly 
Cabibbo-suppressed D+ — > K + ir® decay mode, and a 
new measurement of the D+ — > ir + ir° branching fraction. 

We can compare our results to theoretical expectations 
and evaluate the size of SU(3) violation in these decays. 
In the limit of SU(3), the ratio 



Rsu(3) — 2 



Vcs_ 

Vcd 



T{D+ -» 7T+7T ) 

T(D+ -> K°ir+) 



is expected to approach unity 0. The extra factor of 2 
arises because of the normalization of the 7r° wavefunc- 
tion. Combining values for the CKM matrix elements 
and B{D+ -> K°ir+) taken from @ with our result 



for ZU 



wc find R 



SU(3) 



1.54 ± 0.27. At the 



decay amplitude level this corresponds to a 25% devia- 
tion from S U(3) symmetry, consistent with theoretical 
expectations \U\ . 

We can also compare doubly Cabibbo-suppressed de- 
cays of charged and neutral D mesons. The two decays 
D + — > K + tt° and D° — > K + tt~ differ only in the flavor 
of the spectator quark in the D meson. In the absence 
of D° mixing and taking into account that the D + decay 
includes a 7r° in the final state, the ratio of decay rates is 
expected to be 1/2. This ratio could be modified by W- 
annihilation and W-exchange amplitudes that contribute 
differently to D+ -> K+n and D° -> K+tt~ decays 
Experimentally we find 



T(D+ -» K+tt°) _ B(D+ -> K+TT°) T D 
T(D° -> K+-K-) ~ B(D° -> K+n') t d + 
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